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ABSTRACT 

During its first data cycle, between 2005 and the beginning of year 2006, 
the fast repositioning system of the MAGIC Telescope allowed the observation 
of nine different GRBs as possible sources of Very High Energy (VHE) 7-rays. 
These observations were triggered by alerts from Swift, HETE-II, and Integral; 
they started as fast as possible after the alerts and lasted for several minutes, 
with an energy threshold varying between 80 and 200 GeV, depending upon the 
zenith angle of the burst. No evidence for gamma signals was found, and upper 
limits for the flux were derived for all events using the standard analysis chain of 
MAGIC. For the bursts with measured redshift, the upper limits are compatible 
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with a power law extrapolation, when the intrinsic fluxes are evaluated taking 
into account the attenuation due to the scattering in the Metagalactic Radiation 
Field (MRF). 

Subject headings: gamma rays: bursts — gamma rays: observations 



Introduction 



The physical origin of the enigmatic Gam ma-Ray Bursts is still under debate today, 40 
years after their discovery (see iMeszaros 20061 . for a recent review). The possible detection 
of radiation in the Very High Energy region (extending between few tens of GeV and few 
tens of TeV) will lead to a deeper understanding of the acceleration mechanisms and the 
emission processes from gamma-ray bursts. The 7-ray emission observed by the Energetic 
Gamma-Ray Experiment Telescope (EGRET) in som e case extends up to the VHE band 
( IHurley et al. 1994 iDingus 19951 : iGonzalez et al. 20031 ) , favouring the hypothes is of a highl y 
relativistic source of non-thermal radiation situated in an optically thin region (IPiran 19991 ): 
more insight, however, can be gained by a clear signal detection in the VHE region, or the 
evaluation of stringent upper limit in this energy band. 



Severa l observations of GRBs at energies above 100 GeV have been attempted (IGotting et al. 2003 
Zhou 20031 ). without showing any indication of a signal. This is due to relatively low sensitiv- 
ity, as in satellite-borne detectors, or to high energy thresholds, as in the previous generation 
of Cherenkov telescopes or in particle detector arrays. Only few tentativ e detections of radi- 
ation above 0.1 TeV were reported by MIL AGRITO for GRB 970417a jAtkins et al. 200oh 
and by the GRAND array on GRB 971110 JPoirer et al. 2003h . 



Upper limits on the prompt or delayed emission of GRBs were also set by Whipple 



flConnauehton et al. 19971: iHoran et al. 2007t). MILAGR O (lAtkins et al. 20051: ISaz Parkinson et al. 2006a : 
Saz Parkinson et al. 200 6bl: ISaz Parkinson et al. 2006ch . STACEE djarvis et al. 2005h and 



HEGRA-Airobicc ( 



Padilla et al. 19981 ). 



Imaging Atmospheric Cherenkov Telescopes (IACT) of the latest generation achieve a 
better flux sensitivity and a lower energy threshold, and thus are better suited to detect 
VHE 7-rays; on the other hand, their small fields of view permit unguided observations only 
by virtue of serendipitous detection, and they have to rely on an external trigger, such as 
that provided by automated satellite link to the GRB Coordinates Network (GCN), which 
broadcasts the coordinates of events triggered and selected by dedicated satellite detectors. 



Among the new Cherenkov telescopes, MAGIC (IMirzoyan 20051 ) is best suited for the 
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detection of the prompt emission of GRBs, because of its low energy thresh old, large effective 
area and, in particular, its capability for fast slewing (IBretz et al. 20031 ). The low trigger 
threshold, currently 50 GeV at zenith, should allow the observation of GRBs even at large 
redshift, as lower energy radiation can effectively reach Earth without interacting much with 
the MRF. Moreover, in its fast-slewing mode, MAGIC can be repositioned within 30 s to any 
position on the sky; in case of a target-of-opportunity alert by GCN, an automated procedure 
takes only few seconds to terminate any pending observation, validate the incoming signal, 
and start slewing toward the GRB position. Up to now, the current maximal repositioning 
time is ~ 100 s. In t wo cases, this allowed to put upper limits on the GRB flux e ven during 
the prompt emission (IGalante et al. 20051 ; I Albert et al. 20061 ; iMorris et al. 2006bl ). 



The detection of VHE radiation from the GRB is important for comparing differ- 
ent theoretical models. The emission in the GeV-TeV range in the prompt and delayed 



phase is predicted by several authors (see IZhang fc Meszaros 200 It iPe'er fc Waxman 2004 



Razzaque et al. 2004! . for a detailed analysis). Possible processes com prise leptonic and 

hadronic models: i nverse- Compton (IC) scattering by electrons in internal (IPapathanassoiu fc Meszaros 20C 
Pilla &; Loeb 19981) o r external shocks ( Meszaros et al. 19941) . IC i n the afterglow sh o cks 
rtDermer et al. 200ol : Izhang fc Meszaros 200ll: berishev _etal. 200ll ; IWang et al. 200lh . IC 
by electrons responsible fo r optical flashes ( Beloborodov 2005T). and pur e electron-synchrotron 
(IZhang &; Meszaros 2 00 if); proton-synchrotron emission (jTotani 20001). photon-pion produc- 



tion (IWaxman 19951 ; iBottcher fc Dermer 19981; IChiang fc Dermer 199 



and neutron cascades (IBahcall &; Meszaros 2000 



9 ; 



Derishev et al. 1999 



Li et al. 2002: Fragile et al. 2004T) . 



Rossi et al. 20061 ). 



During the early afte rglow phase, the rec ent observations by the Swift satellite of X-ray 
flares lasting 10 3 4- 10 5 s ( iBurrows et al. 20051 ) suggested an exten ded activity in the cen tral 
engine of the GRB, and thus emission from late internal sho cks (IKobavashi et al. 20071) or 
from refreshed shocks due to energy injections at later time ( iGuetta et al. 20061 ). In some 
cases, the energy release of these flares can be of the same order of magnitude of the energy 
release in the prompt phase, as reported for GRB 050502b. Possibility of correlated 7-ray 
emission extending into the GeV-TeV region is predicted as well, where the correspond- 
i ng VHE flares are predicted to originate from IC scattered photons in the forward shock 
(jWang et al. 20061 ). Thus, observation of the delayed activity is of particular interest being 
in most cases not constrained by the alerting and slewing time, and being still connected to 
the investigation of the central engine dynamics. 

Measurements in this energy range can be used to test all these competing models. 
However, as most of the observed GRBs occur at large redshift, strong attenuation of the 
VHE 7-ray flux is expected, as a result of the interaction with low energy photons of the 
MRF (INikishov 196ll ; lde Jager fc Stecker 20021 ). The knowledge of the redshift, therefore, is 
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important for a precise interpretation ([Mannheim et al. 19961 ). 

In this article, we report on the analysis of data collected on several GRBs followed by 
MAGIC during their prompt emission and early afterglow phases. 



2. Gamma-ray analysis with the MAGIC telescope 



The Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescope (IMirzoyan 20051 ). 
located on the Canary Island of La Palma (2200 m asl, 28°45' N, 17°54' W), completed its 
commissioning phase in early fall 2004. MAGIC is currently the largest IACT, with a 17 m 
diameter tessellated reflector dish consisting of 964 0.5 x 0.5 m 2 diamond-milled aluminium 
mirrors. In its current configuration, the MAGIC photo-multiplier camera has a trigger 
region of 2.0° diameter (ICortina et al. 20051 ). and a trigger collection area for 7-rays of the 



order of 10 5 m 2 , which increases further with the zenith angle of observation. Presently, 
the accessible trigger energy range spans from 50-60 GeV (at small zenith angles) to tens of 
TeV. The MAGIC telescope is focused to 10 km distance - the most likely height at which 
a 50 GeV 7-ray shower has its maximum. The accuracy in reconstructing the direction of 
incoming 7-rays, the point spread function (PSF), is about 0.1°, slightly depending on the 
analysis. 



The reconstructed signals are calibrated (IGaug et al. 20051 ). and then cleaned of spuri- 
ous backgrounds from the light of the night sky using two different image cleaning proce- 
dures: one algorithm requiring signal exceeding fixed reference levels, and a second alg orithm 
employing additionally the reconstructed information of the arrival time (IGaug 20061 ). Non- 
physical background images are eliminated (e.g. car flashes having trig gered the readout) . 
Events are processed by means of t he MAGIC s t andard analysis so ftware (IBretz et al. 20051 ). 
using the standard Hillas analysis (IHillas 19851 ; iFegan et al. 19971 ). Gamma /hadr on separa- 
tion is performed by means of the Random Forest (RF) method (IBreiman 20011 ). a classi- 
fication method that combines sev eral parameters des cribing the shape of the image into a 
new parameter called hadronness ( Hengstebeck 20061 ). the final 7/hadron discriminator in 
our analysis. Monte Carlo samples are used to optimize, as a function of energy, the cuts 
in hadronness. The energy of the 7-ray is also estimated using a RF approach, yielding a 
resolution of ~ 30% at 200 GeV. The parameter alpha of the Hillas analysis, which is related 
to the direction of the incoming shower, is not included in the calculation of hadronness, 
as it is used separately to evaluate the significance of a signal. If the telescope is directed 
at a point-like 7-ray source, as a GRB is expected to be, the alpha-distribution of collected 
photons should peak at 0°, while it is uniform for isotropic background showers. 
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Blind Test with Crab Nebula 



On 2005 October 11 at 02:17:37 UT, the Integral satellite announced GRB 0510110 at 
the position J2000 R.A. = 5 h 34 m 47 s , decl . = +21° 54' 39". A few hours later, Integral sent 



a new GCN flMereghetti fc Mowlavi 20051 ) stating that GRB 051011 was in fact the Crab 
Nebula. Thus, in a blind test, we acquired 2814 seconds of events coming from the Crab 
Nebula, the sta ndard source of 7-r ays at VHE energies. The analysis yielded a 14cx signal 
above 350 GeV (IScapin et al. 20061 ) . showing that MAGIC can observe, at 5a level, spectra 
of 5 Crab Units (1 C.U. = 6.57 x 10 
above 300 GeV, and in 90 s if below 300 GeV. 



10 cm 2 s 1 above 100 GeV) of intensity in 40 s, if 



4. GRBs observed by MAGIC during its first observation cycle 

An automatic alert system has been operational from July 15 th , 2004. Since then, about 
200 GRBs were detected by HETE-II, Integral and Swift, out of which about 100 contained 
GRB coordinates. Time delays to the onset of the burst were of the order of several seconds 
to tens of minutes. During the first MAGIC data cycle, between April 2005 and March 2006, 
9 GRBs were observed by MAGIC during the prompt or the early afterglow emission phase, 
as listed in Table [TJ In two cases the prompt alerting by the GCN and the fast reaction of 
the MAGIC telescope allowed to take data not only on the early afterglow, but also on part 
of the prompt emission of the burst. These two bursts are GRB 0507013a and GRB 050904 
and will be considered separately 



4.1. The properties of observed GRBs 

Table [2] summarizes the properties of observed GRBs by MAGIC according to the in- 
formations distributed through the GCN Circular Service. 



GRB 050421 was detected by the Burst Alert Telescope (BAT) on board Swift (IBarbier et al. 2005al ) 
The other telescope on board Swift, the X-ray Telescope (XRT), could observe the burst in 
the 0.2-10 keV range since T p + 97 s and could detect two X-ray flares at T + 110 s and 



T + 154 s (IGodet et al. 20061 ). Figure [T] shows the X-ray light curve of this burst. It can be 
seen that the MAGIC observation window is overlapped with the XRT ones on the X-ray 
afterglow. In particular, the two small X-ray flares are in the observation window of MAGIC. 



HBAS alert nr. 2673, see at 



http : / / ibas .mi . iasf . cnr . it/ 
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Burst 


Satellite 


T [UT] 


AT a i er t 


AT s tart 


^■slewing 


Data 


zenith angle 


1. 


GRB 050421 


Swift 


04:11:52 


58 S 


108 S 


26 s 


75 min 


~ 52° 


2. 


GRB 050502a 


Integral 


02:13:57 


39 s 


689 s 


223 s 


87 min 


~ 30° 


3. 


c \ T~> t ^ ncAcnc 

GRB 050505 


Swiit 


26:22:21 


540 s 


717 s 


90 s 


101 mm 


A no 

~ 49 


4. 


GRB 050509a 


Swift 


01:46:29 


16 s 


131 s 


108 s 


119 min 


~ 58° 


5. 


GRB 050713a 


Swift 


04:29:02 


13 s 


40 s 


17 s 


37 min 


~ 49° 


6. 


GRB 050904 


Swift 


01:51:44 


82 s 


145 s 


54 s 


147 min 


~ 24° 


7. 


GRB 060121 


HETE-II 


22:24:54 


15 s 


583 s 




53 min 


~ 48° 


8. 


GRB 060203 


Swift 


23:55:35 


171 s 


268 s 


84 s 


43 min 


~ 44° 


9. 


GRB 060206 


Swift 


04:46:53 


16 s 


59 s 


35 s 


49 min 


~ 13° 



Table 1: Summary of GRBs observed by MAGIC from April 2005 to March 2006. AT a i ert 
stands here for the time delay after To until the burst coordinates were received from the 
GCN. ATgtart is the total time delay before the observation could be started, of which Viewing 
is the time lost for repositioning the telescope. Data column shows the total amount of data 
taken. 



No optical counterpart was observed, thus GRB 050421 has been catalogued as a dark burst. 



GRB 050502a was triggered b y Integral, no X-ray counterpart wa s observed, but an 
optical afterglow followed the burst flGotz 2005al botz 2005bl : burig 2005^ . GRB 050505 was 
triggered by BAT and its light curve presented three s hort spikes at T +23.3 s, T +30.4 s and 
T +50.4s (IHurkett et al. 2005al : iHullinger et al. 20051 ). Both X-ray and optical observations 
followed the burst, but there was no simultaneous observation by MAGIC and the other 
instruments on board the satellite, as shown in figure [2j 

GRB 050509a was triggered by BAT a nd later XRT could detect an X-ray counterpart 
Jrlurkett et al. 2005bl : barbier et al. 2005bh . GRB 060121 is the only short burst observe d 
by MAGIC and was triggered by HETE-II jArimoto et al. 20061 : iGolenetskii et al. 20061 ) . 
XRT obser ved the afterglow as w ell and detected a fading X-ray source inside the HETE II 
error box (IMangano et al. 20061 ). as shown in figure [31 Unfortunately, there is no over- 
lap between MAGIC observation and XRT observation. An optical counterpart was not 
confirmed by TNG, which did detect, however, a weak source inside the XRT error box 
(IMalesani et al. 20061 ). Moreover, HST gave no evidenc e of an optical aft erglow, although 
the burst lay close to a faint red galaxy at high redshift (ILevan et al. 20061 ) . 



T he last two bursts, GRB 060203 and GRB 060206, are both long bursts triggered by 
BAT faarthelmv et al. 20061 : iMorris et al. 2006ah . MAGIC data overlap with XRT data on 
the X-ray afterglow of GRB 060206, immediately after BAT data as shown in figure HI No 
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Burst 


Trigger # 


Energy Range 


^90 


Fluence 


Z 


1. 


GRB 


050421 


115135 


15-350 keV 


10 s 


1.8 x 


10" 


7 


- 


2. 


GRB 


050502a 


2484 


20-200 keV 


20 s 


1.4 x 


10" 


-6 


3.79 


Q 
O. 


GRB 


050505 


11 / t)U4 


lo-oou Rev 


ou s 


4.1 x 


10" 


6 


4.Z I 


4. 


GRB 


050509a 


118707 


15-350 keV 


13 s 


4.6 x 


10" 


7 




5. 


GRB 


050713a 


145675 


15-350 keV 


70 s 


9.1 x 


10" 


6 




6. 


GRB 


050904 


153514 


15-350 keV 


225 s 


5.4 x 


10" 


6 


6.29 


7. 


GRB 


060121 


4010 


0.02-1 MeV 


2 s 


4.7 x 


io- 


-6 




8. 


GRB 


060203 


180151 


15-350 keV 


60 s 


8.5 x 


io- 


7 




9. 


GRB 


060206 


180455 


15-350 keV 


11 s 


8.4 x 


io- 


7 


4.05 



Table 2: Main properties of GRBs observed by MAGIC. The third column shows the typical 
energy range of detector on board of the satellite, while in the fourth, fifth and sixth column 
there are the corresponding measured duration Tgo, fluence in [erg cm -2 ] and redshift. 



evidence of flares or of the jet break have been claimed, but optical observations provided a 
high redshift value (jPalmer et al. 20061 : lAoki et al. 20061 ) . 



4.2. GRB 050713a and GRB 050904, prompt emission observation 



GRB 050713a is of particular interest, being the first burst observed by MAGIC dur- 
ing its prompt emission jAlbert et al. 2006h . On July 13 th , 2005 at 04:29:02 UT the BAT 
instrument detected a burst located at R.A. = 21 h 22 m 09!53, decl. = +77° 04' 29"50 ± 3' 



( jFalcone et al. 20051 ). The MAGIC alert system received and validate d the alert 13 s afte r 
the burst, data taking started 40 s after the burst original time T ( IGalante et al. 20051 ). 
The burst was classified as a bright burst by Swift with a duration of T90 = 70 ± 10 s. The 
brightest part of the keV emission occured within T + 20 s, three smaller peaks followed 
at T + 50 s, T + 65 s and T + 105 s, while a pre-burst peak took place at T — 60 s (see 
figure EJ). The spectrum, over the interval from T — 70 s to T + 121 s, can be fitted with a 
power law with ph oton index —1.58 ± 0.07 and yields a fluence of 9.1 x 10~ 6 erg cm -2 in the 
15-350 keV range (IPalmer et al. 20051 ). The burst triggered also the Konus-Wind satellite 
(IGolenetskii et al. 20051 ). which measured the spectrum of the burst during the first 16 s, 
that is the duration of the first big peak as reported by Swift. In the local coordinate system 
of MAGIC, GRB 050713a was located at an azimuth angle of —6° (near North) and a zenith 
angle of 50°. The sky region of the burst was observed during 37min, until twilight. 



Also GRB 050904 is of particular interest, being the second and the latest burst with 
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Fig. 1. — Flux of GRB 050421 measured by BAT and XRT. The orange shadowed area 
represents the MAGIC observation time window and the overlap with Swift data. 



prompt emission observed by MAGIC. It was triggered at 01:51:44 UT by BAT, coordi- 
nates were R.A. = h 54 m 50 s . 79, decl. = +14° 05' 09".42 ± 3' jCummings et al. 2005h . XRT 
slewed promptly and started the observation at T + 161 s, revealing an uncatalogued fading 
source. It is a lon g burst (T 3 n = 225 s), with a total fluence of 5.4 x 10~ 6 erg cm -2 in the 



15-150 keV range (jSakamoto et al. 20051). This burst is the most distant burst ever observed, 
with an estimated redshift z = 6.29 (IKawai et al. 20051). I ts X-ray light curve (see figure [6]) 
shows a clear X-ray flare at T + 466 s (jMineo et al. 20051 ). thus in the MAGIC observation 
window. 



5. Results 

All nine GRBs were analyzed using the MAGIC standard analysis described above. In 
this work the image cleaning algorithm using also arrival time information was used, as it 
is more robust. For each GRB a dedicated OFF-source data set was selected on the basis 
of being compatible with the ON-data with respect to several parameters: zenith angle, as 
the effective area depends strongly on it; local brightness of the sky, depending mostly on 
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Fig. 2. — Flux of GRB 050505 measured by BAT and XRT. The orange shadowed area 
represents the MAGIC observation time window, starting 717 s after the burst onset. 



Moon phase and zenith angle; trigger rate, depending mostly on atmospheric transparency 
and on hardware settings. Loose preliminary cuts were used to remove unphysical events. 
After training of the RF, for each burst a hadronness cut was applied which could grant 
about 90% efficiency on 7-ray events according to the corresponding Monte Carlo, in order 
to keep high statistics of possible 7-rays events. 

The analysis showed no evident signal excess, as can be seen in figure [7J For each burst 
the alpha plot over the whole data set and for reconstructed energies greater than 100 GeV is 
shown. The alpha distributions of the GRB data sets are flat, as expected from background 
hadronic events, and are compatible with the corresponding OFF-source data set. No excess 
in the signal region, i.e. for alpha < 30°, can be seen. 

No excess was also evidenced using a temporal analysis: The entire data taking interval 
was divided into 20-second time bins and the number of potential 7-ray events, extracted 
from the alpha distribution; they are shown in the light curves of Figure [SJ red filled circles 
denote the excess events, blue open circles the background events (offset by 5 from excess 
counts in order to make the plot more readable). The distributions of excess events remain 
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Fig. 3.— Flux of the afterglow of GRB 060121 measured by Swift XRT detector. MAGIC 
observation window is shown as orange area. 



zero on average during the observation, and no significant variation of the sample average is 
visible with time. 

Upper limits have been derived for the first 30 minutes of each burst using the Rolke 
approach (IRolke et al. 20051 ) in six reconstructed energy bins: 80-120 GeV, 120-175 GeV, 
175-225 GeV, 225-300 GeV, 300-400 GeV and 400-1000 GeV. A systematic uncertainty of 
30% on the efficiency has been considered in the upper limit calculation. In every recon- 
structed energy bin the upper limit in number of excess, calculated with the Rolke approach, 
has been converted in flux unit using the effective 3X6 cl, clS explained in appendix |A] The 
typical effective area of MAGIC for different zenith angle is shown in figure M Table [3] 
summarizes the upper limits derived for all nine GRBs during the first 30 minutes of data 
taking. 



6. Discussion 

A preliminary estimate of the observability of GRB s by the MAGIC tel e scope had orig - 



inally been derived using the fourth BATSE catalogue (IBastieri et al. 20051 ; iGalante 20061 ). 
The GRB spectra were extended to GeV energies with a simple power law and using the 
observed high-energy spectral index. The extrapolated fluxes were finally compared to the 
estimated MAGIC sensitivity. Setting conservative cuts on observation times and signifi- 
cances, and assuming an energy threshold of 30 GeV, a 5 cr-signal rate of 0.5 — 2 per year was 
obtained for an assumed observation delay between 15 and 60 s and a BATSE trigger rate of 
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Energy bin 
[GeV] 


Energy 
[GeV] 


Fluence Upper Limit 
[cm 2 keV -1 ] [erg era" 2 ] C.U. 


GRB 060121 


120-175 
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2.64 x HT 15 9.67 x 10" 8 0.41 


175-225 


212.8 


6.57 x 10~ 16 4.76 x 10" 8 0.25 


225-300 


273.7 


2.13 x 10~ 16 2.56 x 10" 8 0.15 


300-400 


367.7 


4.47 x 10~ 16 9.66 x 10" 8 0.69 


400-1000 


636.4 


4.84 x 10~ 17 3.14 x 10~ 8 0.31 


GRB 060203 


120-175 


151.5 


1.10 x 10~ 14 4.03 x 10~ 7 1.71 


175-225 


219.5 


5.07 x 10~ 16 3.91 x 10" 8 0.21 


225-300 


274.0 
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300-400 
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400-1000 
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GRB 060206 


80-120 


85.5 


1.23 x 10~ 14 1.44 x HT 7 0.44 


120-175 


139.9 


9.83 x 10~ 16 3.08 x 10" 8 0.13 


175-225 


210.3 


5.50 x 10~ 16 3.89 x 10" 8 0.20 


225-300 


269.2 


3.65 x 10~ 16 4.23 x 10" 8 0.25 


300-400 


355.4 


6.47 x 10~ 16 1.31 x 10~ 7 0.91 


400-1000 


614.0 


2.88 x 10~ 17 1.74 x 10~ 8 0.17 



Table 3: Derived fluence upper limits for the first 30 minutes of data of nine Gamma- Ray 
Bursts. The first column shows the reconstructed energy bins in which the analysis has 
been done. The second column shows the true energy at which the upper limits have been 
calculated, and is the energy giving the average flux upper limit in the reconstructed energy 
bin. The last column shows the upper limit value in Crab Unit. 
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Fig. 4. — Light curve of GRB 060206 measured by BAT, sampled in bins of 1 s. The 
beginning of MAGIC observation is shown shadowed. 



'360/year. Taking into account the rate of GRBs (IGuetta et al. 20031 ). and extrapolating 



GRB spectra , as observed by B ATSE, to VHE using an unbroken power law of reasonable 



power index (IPreece et al. 2000l ). it was foreseen that MAGIC could detect about one GRB 



per year at a 5<r level. A maximal redshift up to z = 2 was considered. 

This estimate must be revised: the Swift alert rate is about a factor 2 lower than 
predicted, and it includes bursts even fainter than those observed by BATSE; also, for these 
bursts the effective MAGIC energy threshold at analysis level was higher than the assumed 
one (~ 80 GeV); most important, the distribution of bursts detected by Swift has a much 
higher median redshift than expected. As a result, the number of GRBs that MAGIC can 
detect is now estimated to lie in the range of 0.2-0.7 per year. This number can be expected 
to increase again with the launch of the GLAST and AGILE satellites, and the increased 
number of alerts due to the GRB monitoring by GLAST, AGILE and Swift altogether. 
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Fig. 5. — Flux of the prompt and afterglow emission of GRB 050713a measured by BAT 
and XRT. The orange shadowed area illustrates the observation time window of MAGIC. 



7. Conclusions 



MAGIC was able to observe part of the prompt and the early afterglow emission phase of 
many GRBs as a response to the alert system provided by several satellites. No excess events 
above ~ 100 GeV were detected, neither during the prompt emission phase nor during the 
early afterglow. We have derived upper limits for the 7-ray flux between 85 and 1000 GeV. 
These limits are compatible with a naive extension of the power law spectrum, when the 
redshift is known, up to hundreds of GeV. 

For the first time an Atmospheric Cherenkov telescope was able to perform direct rapid 
observations of the prompt emission phase of GRBs. This is particulary relevant in the 
so called "Swift era". Although strong absorption of the high-energy 7-ray flux by the 
MRF is expected at high redshifts, given its sensitivity to low fluxes and its fast slewing 
capabilities, the MAGIC telescope is currently expected to detect about 0.5 GRBs per year, 
if the GRB spectra extend to the energy domain of hundreds of GeV, following a power law 
with reasonable indices. 
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Fig. 6. — Flux measured by BAT and XRT in the prompt and afterglow emission of 
GRB 050904. The MAGIC observation window is shown shadowed. 
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A. Upper Limit Calculation 

The recorded number of events in a particular reconstructed energy bin AE vec is: 

POO 

AW = / <p(E)A cS (E\AE vcc ) dE x AT (Al) 
Jo 




Fig. 7. — Alpha plots of all nine GRBs for the complete data set of each burst. Points with 
error bars refer to the burst data, the line refers to the background. 
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Fig. 9. — Effective area of MAGIC, after typical cuts used in this analysis, for three different 
zenith angles. 
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where <p(E) is the flux (ph cm -2 s -1 GeV -1 ), A e g(E\ AE rec ) is the effective area after all cuts, 
included the reconstructed energy cut AE TCC , and AT the total time interval of observation. 
It should be noted that the flux <p and the effective area A depend on the true energy, while 
the cuts for the selection of the excess events N AEicc and of the effective area A eS (E\AE rec ) 
depend on the reconstructed (estimated) energy. The integral is computed in true energy 
dE. 

Being the effective area A eS (E\ AE rec ) depending on energy, we must assume a spectral 
shape, in our case the typical power law of a GRB: 



m = k- (j^ (A2) 

where k is a normalization factor (ph cm -2 s -1 GeV -1 ), (3 is the average high energy power- 
law index (3 = —2.5 and Eq the pivot energy (e.g. 1 GeV). From equation ( lAlf) we obtain 
the normalization factor as 



* A e g(E\AE Tec )(E/E o y dE x AT 1 

In our case N A e icc is the upper limit in number of excess events calculated with Rolke 
statistics. The flux upper limit is 



UL (T) = K ■ I- (A4) 

The value of the energy E is chosen in a very simple way. We calculate the energy which 
$ives the average flux in the observed reconstructed energy bin: 



f~0(T)A cff (T|AT rcc )dT _ N AEicc 
WAcff / °° A cH (E\AE rec ) dE J °° A cS (E\AE rcc ) dE x AT 1 ° j 

From equations (lAlj) and (1A3j) we can write 



N AErcc J™A eS (E\AE Tec )(E/E o ydE 

r °° A eS (E\ AE TCC )(E/E )P dE ' A eS (E\AE Tec ) dE x AT 



K ■ ((E/Eof)^ 



(A6) 
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Denning (E*/E ) = {(E/E^f) 1 ^ from equation flAil) 



we can calculate the average flux 



upper limit in the reconstructed energy bin: 



</>ul(£*) = h 




cff 



(A7) 



Equation IA7I has been used to calculate the upper limits shown in table [3j 
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